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METHODS  AND  PROGRAMS  TO  UTILIZE  SATELLITE  INDIRECT  SOUNDINGS 
IN  ENVIRONMENTAL  ANALYSIS  AND  PREDICTION 


1.  INTRODUCTION 

1 . 1  The  Assignment 

The  general  objectives  are  methods  and  programs  for  utilizing  remote 
soundings  of  the  atmosphere  from  space  fcr  improving  environmental  analysis 
and  prediction.  This  technical  report  describes  the  development  of  a 
capability  for  diagnosing  radiances  measured  by  a  multi-channel  scanning 
radiometer,  for  inherent  clear-column  radiance  estimates,  and  a  method 
for  using  this  diagnosed  information  for  enhanced  resolution  of  the 
atmospheric  thermal- structure  distribution.  Conceptions  and  formulations 
are  described;  complete  details  appear  elsewhere,  with  the  program 
documentation. 

This  project  is  a  continuation  of  earlier  work  also  performed  by 
Meteorology  International  Incorporated.  The  pertinent  concepts  and 
formulations  for  the  diagnosis  of  measured  radiances,  for  the  inherent 
resolution  of  clear-column  radiances,  first  appeared  in  a  processing  scheme 
which  we  denoted  by  the  letters  CLRX.  The  scheme  is  based  on  an 
idealization  of  cloud  structures  and  hence  may  be  termed  a  diagnostic  model. 
The  CLRX  model  was  designed  for  the  diagnosis  of  individual  scanlines 
measured  by  Special  Sensor  E  (SSE)  aboard  Defense  Meteorological  Satellite 
Program  (DMSP)  spacecraft.  A  later  version  of  the  model,  named  CLRXN, 
was  developed  for  application  to  the  Vertical  Temperature  Profile  Radiometer 
(VTPR)  sensor  aboard  NOAA  spacecraft,  again  designed  for  diagnosing  single 
scanlines  but  with  exploitation  of  the  greater  proximity  of  adjoining 
scanlines . 
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The  SSE  is  a  radiometer  designed  to  measure  thermally-emitted 


infrared  radiation  from  the  earth-atmosphere  system.  Figure  1  is  an  example 
^  of  the  scan  pattern  of  the  SSE.  The  instrument  scans  from  left  to  right  in 

steps,  providing  25  observation  views  (i.e.  ,  spots)  along  each  scanline 
crosswise  to  the  motion  of  the  satellite.  The  surface  field  of  view  is  about 
20  nm  in  diameter  directly  below  the  satellite,  becoming  broadened  and 
elongated  at  either  end  of  the  scanline.  Eight  radiation  measurements  are 
made  at  each  scan  spot.  SLx  of  the  eight  radiation  observations  are 
provided  by  SSE  channels  positioned  at  different  wavelength  intervals  in 
.  the  15-u“  absorption  band  of  carbon  dioxide.  The  seventh  SSE  channel  is 

located  in  the  12-u-  atmospheric  window  and  the  eighth  channel  senses 
radiation  in  a  water-vapor  absorption  band  near  18  un  . 

The  VTFR  instrument  is  similar  in  complement  of  channels  but  differs 
in  scanning  geometry.  The  VTFR  system  and  processing  job  stream  have 
been  described  by  McMillan  et  al.  (1973). 

Specific  assignments  of  this  project  are 

(i)  to  improve  and  refine  the  diagnostic  model  to  extend  its 
applicability  and  to  increase  the  information  yield,  and 

lii)  to  organize  the  model  into  a  programmed  job  stream  which 
is  capable  of  being  run  routinely  on  the  computer  system  of 
the  Fleet  Numerical  Weather  Central  (FNWC)  for  exploiting 
the  SSE  sensor  of  the  PMS?  satellite  system. 

ft 

These  assignments  have  now  been  performed.  When  implemented 
into  operations  the  system  should  have  a  major  impact  in  advancing 
atmospheric  analysis  and  prediction  capabilities ,  provided  that  the  channel 
source  functions  are  properly  calibrated  and  that  the  new  diagnostic  model, 
named  CLR3LCK,  is  adequately  tuned.  This  realization  will  provide  a  basis 
for  further  improvements  in  the  design  of  satellite  sounding  systems, 

9  respecting  sensors  and  scanning  geometries,  and  for  further  operational 

advances  in  information  assimilation  and  blending. 


1.2 


Earlier  Work 


1.2.1  Specific  Developments  Performed  for  the  Naval 

Environmental  Prediction  Research  Facility  (NEPRF) 

Meteorology  International  Incorporated  (Mil)  performed  Phases  Cr.e 
and  Two  of  the  development  of  techniques  for  operational  exploitation  of 
satellite-observed  sounding  data  for  resolution  of  atmospheric  thermal- 
structure  variabilities.  Developments  under  Phase  One,  reported  by  Hell 
(1972,  1973a,  1973b),  included  a  general  scheme  for  blending  dear-column 
radiance  estimates,  weighted1  as  thermal-structure  information,  into 
enhanced  resolution  of  the  thermal  structure.  This  scheme  is  an  example 
of  what  we  call  structure  blending.  The  subject  version,  named  STRBLND , 
also  exploits  the  information  that  the  atmospheric  is  statically  stable  in 
the  object  scale  of  resolution. 


The  major  development  of  Phase  Two,  reported  by  Hell  (1975a),  is 
a  oasic  scheme  named  CLRX  for  processing  a  scanline  of  multi-channel 
views.  The  scheme  estimates  which  measured  radiances  are  effectively 
free  of  clouds,  and  diagnoses  clear-column  radiance  components  in  regions 
of  cloud-contaminated  views.  The  diagnosed  yield  is  limited  to  those 
channels  which,  in  penetrating  down  to  cloud  tops  in  a  region,  also  include 
some  views  which  are  a:  least  partially  clear  to  the  depth  of  channel 
penetration.  This  basic  version  of  CLRX  contains  the  essential  formulations 
with  applicability  limited  to  oceanic  regions  having  little  spatial  variability 
in  sea-surface  temperature.  There  are  additional  limitations  and  a 
complexity  of  contributing  variances. 


In  the  FI3  methodology  developed  by  Mil  all  parameter  estimates 
are  generally  weighted.  The  weight  of  an  estimate  is  equivalent  to  the 
inverse  of  the  associated  variance.  The  weight  may  refer  to  the  reliability, 
or  independent  worth,  of  a  parameter  estimate  in  a  particular  context,  or  to 
the  resolution  accuracy  of  an  estimate  resulting  from  a  blending  of  assembled 
information. 


The  information  inherent  in  measured  radiance  scans  is  highly  variable. 
It  is  related  to  the  characteristics  of  the  cloudiness  that  is  present.  The 
diagnostic  yield  generally  decreases  from  channel  deeper-seeing  channel. 
The  resolution  weights  range  from  minimal  for  channels  fully  blocked  by 
cloud  in  a  region,  to  very  high  in  clear  regions.  Because  of  this  great 
variability  in  information  yield,  the  weight  associated  with  a  diagnosed 
clear-column  radiance  estimate  is  an  essential  component  of  the  CLRX- 
diagnosec  information.  Because  of  this  great  variability  in  resolution,  the 
diagnosed  clear-column  radiances  cannot  be  categorically  forced  on  the 
atmospheric  thermal-structure  analysis.  Furthermore  the  resolution  weights 
for  the  diagnosed  clear-column  radiances  must  be  reduced  for  the  associated 
channel-calibration  variances  if  they  are  to  be  blended  with  other  independent 
estimates  and  expressions  of  information  bearing  directly  on  the  atmospheric 
thermal  structure.  The  structure-blending  scheme  developed  under  Phase 
One  encompasses  these  considerations. 

Subsequent  to  Phase  Two,  Mil  tested  the  basic  version  of  CLRX  on 
a  sample  of  scanlines  measured  by  sensors  of  two  DMSP  satellites.  Despite 
limitations  of  the  available  data  sample  (e.g.,  incomplete  scanlines,  a 
region  of  islands  and  sea-surface-temperature  irregularities,  and 
uncertainties  in  calibration  relative  to  calculated  first-guess  radiance 
estimates)  these  tests  demonstrated  the  functional  capability  of  CLRX. 

The  results  were  reported  by  Holl  (1975b).  The  wide  range  in  yield  was 
quite  apparent. 


1.2.2  Diagnostic  Model  CLRXN 

The  development  of  an  improved  version  of  CLRX,  designed  for  the 
scanning  geometry  of  VTPR  systems  aboard  NOAA  spacecraft,  was  reported  by 
Holl  et  al.  (1976a).  The  VTPR  scanlines  are  not  as  long  in  surface  nadir-ancle 
sweep,  and  are  about  three  times  more  frequent  along  the  orbit  than  are  the 
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SSE  scanlines.  In  applying  :he  model. spots  are  paired  for  contrasting  local 
radiance  variabilities.  In  CLRXN  these  pairings  are  not  limited  to  being  of 
the  same  scanline.  Three  successive  scanlines  are  analyzed  for  information 
elements  used  in  the  diagnosis  of  the  clear-column  radiance  components  of 
the  middle  scan.  Several  improvements  were  introduced  in  CLRXN. 

Following  delivery  of  the  CLRXN  software  package  to  the  National 

Environmental  Satellite  Service  (NESS)  of  NCAA,  that  organization  evaluated 
2 

the  capabilities."  The  data  base  that  was  used  for  the  test  is  a  simulated 
package:  VTPR  radiances  calculated  for  a  variety  of  cloud  structures.  The 
CLRXN  capabilities  were  considered  to  be  superior  to  other  procedures 
compared  in  the  tests.  Seme  advantages  and  some  disadvantages  of  CLRXN 
were  listed.  Four  advantages  were  stated  to  be  distinct:  CLRXN  produces 
an  independent  estimate  of  the  clear-column  component  of  the  window 
channel:  a  weighted  clear-column  estimate  is  produced  for  each  channel 
and  spot,  and  hence  the  potential  information  yield  is  greater;  the  value 
of  the  associated  weights  is  recognized;  the  program  can  be  run  without 
independent  sea-surface-temperature  estimates  and  hence  its  extension  to 
land  areas  is  almost  immediate  for  flat  homogeneous  terrain. 

Three  of  the  four  listed  disadvantages  can  readily  be  discounted  or 
remedied:  The  apparent  dependence  on  independent  (i.e.,  first-guess 
calculated)  radiance  estimates  is  adjustable;  clear-column  estimates  are 
not  produced  for  spots  at  scanline  extremes  (i.e.,  spots  1,  2  and  22,  23) 
because  the  potential  yield  dees  not  warrant  the  programming  details:  a 
ciear-coiumn  estimate  of  the  water  vapor  channel  is  not  produced  but  this  is 
recognized  as  a  simple  adjunct  to  the  program — the  yield  will  of  course  be 


2 

“In-house  memorandum  by  Dr.  Henry  Fleming  of  the  National 
Environmental  Satellite  Service,  NCAA,  Department  of  Commeroe,  dated 
September  16,  19"6.  Information  oepv  sent  to  MIL 
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dependent  on  the  local  cloudiness  characteristics.  The  fourth  disadvantage 
mentioned  is  that  CLRXN  requires  about  7  or  8  times  as  much  data  processing 
as  do  the  schemes  that  NESS  has  been  using  for  this  Component  in  their  job 
stream.  This  processing  cost,  supplemental  to  the  costs  of  all  other  related 
components  of  the  satellite  and  ground  processing  systems,  must  be  weighed 
against  the  additional,  and  potential,  yields  of  information. 

Another  item  that  is  apparent  in  the  NESS  evaluation  is  that  the  errors 
in  the  clear-column  radiance  components  diagnosed  by  CLRXN  tend  to  be 
negative;  that  is,  the  results  favor  lowering  the  correct  radiance  values. 

The  diagnosed  results  are  consistent  with  the  associated  weights,  and 
there  is  no  way  to  interpret  this  error  inclination  as  a  calibration  bias. 

We  now  believe  this  bias  to  be  due  to  the  diagnostic  model  as  applied  in 
both  CLRX  and  CLRXN:  The  cloud  structure  is  interpreted  as  being  composed, 
locally,  of  cloud  having  broken  tops  at  only  one  level  in  the  atmosphere 
and  that  in  the  breaks  the  atmosphere  is  clear.  Although  this  diagnosed 
information  is  weighted  according  to  assessed  model  conformance  the 
diagnosed  adjustment  axis  is  imposed  on  the  whole  set  of  radiances . 
Disparities  due  to  more  complex  cloud  structures  generally  bias  the  results 
in  the  negative  direction. 


All  of  the  above  considerations  have  been  addressed  in  the  present 
project.  The  new  version  of  the  diagnostic  model  processes  a  block  of 
scanlines  in  each  application.  The  name  CLRBLOK  seemed  appropriate. 
The  major  improvements  and  refinements  are  listed  in  Section  4.1. 


1.2.3  The  Development  of  Blending  by  Weighted  Spreading 


The  major  arena  in  which  the  Fields  by  Information  31ending  1F13) 
methodology  was  conceived  and  advanced  has  been  in  our  development  of 
objective  analysis  techniques  for  a  diversity  of  two-dimensional  fields. 
Other  unique  contributions  have  arisen  in  structure-blending  applications 
and  in  the  subject  diagnostic-model  developments.  The  methodology  has 
a  unity  which  applies  to  all  processing  of  real  Information. 

Since  the  development  of  CLRXN,  Mil  has  produced  a  new  analysis 
capability  designed  for  the  horizontal  distribution  of  atmospheric  upper-air 
parameters,  reported  by  Holl  (1976b).  In  the  course  of  this  work  we 
discovered  an  advancement  which  is  of  major  significance  for  the  solution 
of  equations  by  the  method  of  successive  approximations .  The  new  method 
is  specific  to  the  blending  component  of  all  applications  of  the  FIB 
methodology.  We  refer  to  this  new  method  as  Blending  by  Weighted 
Spreading  (3WS) . 

The  BWS  method  is  characterized  by  the  association  of  a  weight 
with  each  successive  approximation  of  each  element  of  the  developing 
solution.  This  weight  is  an  absolute  measure  of  the  proper  ion  of  the  total 
due  influence  that  has  arrived  at  that  element.  It  is  not  a  measure  of  the 
firmness  of  the  solution  at  that  point  in  the  Iterations,  for  further  Influence 
may  subsequently  arrive  to  substantially  alter  the  developing  value. 
However  when  the  Iterations  are  performed  successively,  in  orderings 
which  alternate  direction  in  each  successive  pass,  most  of  the  due 
influence  Is  felt  in  only  a  few  passes. 

The  3WS  method  does  not  make  use  of  any  first -guess  estimate  to 
the  solution  for  the  purpose  of  initiating  the  successive  approximations . 

The  initiation  stems  entirely  from  whatever  absolute  estimates  have  been 
assembled  at  one  or  more  elements.  Hence  there  is  no  bias  introduced 
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into  the  solution;  characteristic  components  of  such  biases  can  linger  even 
after  many  iterations  using  Successive  Over-Relaxation  (SOR)  methods. 

The  BWS  solution  is  entirely  free  of  extraneous  error  at  every  stage  of  the 
iterations . 

Although  weighted  spreading  ultimately  approaches  the  exact  solution 
asymptotically,  the  initial  iterations  are  amazingly  effective  in  gathering 
information,  when  using  successive  alternating-direction  orderings.  The 
new  diagnostic  model,  CLRBLOK,  introduces  a  further  refinement  which 
exploits  the  best  features  of  both  methods,  BWS  and  SOR.  The  iterations 
begin  with  the  BWS  method.  Then,  wherever  and  whenever  the  accruing 
influence  weight  of  an  element  passes  a  specified  threshold,  the  algorithm 
shifts  into  an  SOR  iteration  for  that  element.  As  the  influence  weights 
approach  unity  they  become  a  redundant  burden  in  the  calculations;  they 
are  then  dropped  for  advanced  elements  and  an  over-relaxation  factor  is 
introduced . 

It  is  the  effectiveness  and  speed  of  the  BWS/SOR  scheme  which 
makes  it  feasible  to  expand  the  capabilities  of  the  new  diagnostic  model, 
CLRBLOK,  for  greater  information  yield  in  quality  and  quantity. 


1.3.1  The  Ideal 


The  atmospheric  conditions  that  are  most  amenable  tc  resolution  of 
clear-column  radiance  components — by  diagnosis  of  suitable  radiance 
measurements  from  space — are,  of  course,  local  regions  that  are  cloud  free. 
But,  in  addition  to  being  cloud  free,  the  atmospheric  thermal  structure  should 
have  low  spatial  variance;  this  is  necessary  for  assessing  that  the  region  is 
indeed  relatively  cloud  free. 
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The  most  amenable  of  the  cloud  conditions  begins  with  a 
scattered-to-broken  cloud  layer  having  a  uniform  top  at  one  distinct  level. 
This  level  is  best  positioned  in  the  midst  of  the  several  channel  source 
functions  of  the  satellite  multi-channel  sensor.  The  scale  of  the  cloud 
pattern  (i.e. ,  the  dimensions  of  the  clouds  and  the  clear  breaks)  should  be 
small  compared  to  the  scale  of  the  spatial  variability  of  the  clear-column 
and  the  cloud-top  radiance  components.  We  refer  to  this  situation  as  a 
two-state  mix  of  radiance  components. 

The  next  most  amenable  of  cloud  conditions  is  a  local  region  in 
which  there  are  two  distinct-and-steadv.  scattered-to-broken,  cloud-top 
levels— a  three-state  mix  of  radiance  components. 

Theoretically  the  tractability  extends  to  more-than-two  clcud-tcp 
levels,  but  as  the  cloud  structure  becomes  more  complex  the  required 
conditions,  of  distinct  and  steady  radiance  components,  are  rarely  to  be 
found  in  the  atmosphere. 


The  satellite  sensing  system  which  is  necessary  for  the  resolution 
of  these  ideal  atmospheric  conditions  is  one  which  measures  the  radiances 
in  distinct  channels,  in  each  field  of  view.  The  more  channels  the  better, 
and  the  scanning  geometry  must  produce  a  high  space-and-time  density  of 
multi-channel  views — dense  in  comparison  to  the  spatial  variability  of  the 
individual  radiance  components.  The  channel  source  functions  may  overlap 
in  the  clear  atmospheric  column,  but  they  should  be  different  from  one 
another  to  the  extent  allowable  by  the  number  of  channels.  The  individual 
multi-channel  views  (i.e. ,  spots)  may  partially  overlap  one  another  but 
the  spots  should  be  small  compared  to  the  spatial  variability  of  the  individual 
radiance  components. 


Given  an  adequate  set  of  channels,  denoted  by  V  in  number  of 
channels,  and  an  adequate  sample  of  views  in  a  small  region,  then  the 
specific  mathematical  procedure  for  analyzing  the  local  cloud  structure  is 
Factor  Analysis.  Many  meteorologists  are  familiar  with  Lorenz's  use  (1956) 
of  Factor  Analysis  to  define  his  Empirical  Orthogonal  Functions.  [Psychologist 
in  particular,  are  mere  familiar  with  the  Factor  Analysis  terminology.]  Factor 
Analysis  has  also  spread  into  many  other  contexts  in  Meteorology. 

Factor  Analysis  assesses  the  variability  of  the  sample  in  terms  of 
orthogonal  axes  in  V  space,  and  determines  the  sample  variance  which  is 
explained  by  each  axis  or  combination  of  axes.  A  sample  with  low  variance 
implies  a  cloud-free  structure  or  a  clear  structure  above  a  uniform  solid 
overcast.  If  one  axis  (i.e. ,  factor)  effectively  explains  most  of  the  sample 
variance  then  the  interpretation  is  two  state:  Each  vector  of  the  sample 
(i.e. ,  each  spot  of  V  radiances)  lies  close  to  a  line  in  V  space.  The  two 
radiance  components — clear  and  (presumed)  cloud  free — bracket  all  the 
measured  mixes.  If  two  axis  (i.e. ,  factors)  effectively  explain  most  of 
the  sample  variance  then  the  interpretation  is  three  state:  All  vectors  of  the 
sample  lie  close  to  a  plane  defined  by  the  two  axes  in  V  space.  The 
orientations  of  the  two  axes  within  this  plane  are  irrelevant.  The  three 
radiance  components  define  the  corners  of  the  triangle  in  which  all  mixes 
fall. 

In  interpreting  a  two-state  cr  a  three-state  mix,  the  residual  variances 
are  a  measure  of  the  difference  between  the  actual  cloud  structure  and  the 
assumed  diagnostic  state.  The  weights  to  be  placed  on  components  of  the 
adjustment  axis  or  adjustment  plane  are  related  to  these  residual  variances. 


The  diagnostic  model  basically  consists  of  diagnosing  inherent 
elements  of  Information  for  each  spot,  and  of  blending  these  elements  for 
the  resultant  determination  of  the  dear-column  radiance  components  for  all 
channels  of  the  continuum  of  spots. 

The  information  produced  by  Tactcr  Analysis  includes  the  assessment 
that  some  measured  radiances  may  be  cloud  free.  This  is  based  on  the 
determination  of  very  low  variance  for  channels  in  the  ambience  of  the  spot. 
The  subset  of  channels  which  arc  assessed  as  possibly  being  cloud  free 
for  any  single  spot  must  be  mutually  consistent  with  the  distribution  of  their 
atmospheric  source  functions.  This  consistency  requirement  should  be 
exploited  in  the  diagnosis  of  the  information  elements. 

The  condition  which  enables  the  blending  of  all  diagnosed  information 
elements  into  a  continuum  Is  that  the  atmospheric  thermal-structure  has  a 
high  degree  of  horizontal  continuity  in  the  object  scale  of  resolution.  The 
diagnostic  model  imposes  a  degree  of  continuity  on  the  dear-column 
atmospheric  radiance  components  from  spot  to  spot.  The  degree  of  validity 
of  this  continuity  imposition  can  be  enhanced  by  pre-processing  the  measured 
radiances  before  diagnosis.  These  preprocessing  components  also  enhance 
the  cloud-structure  diagnosis. 


The  efficacy  of  remote  sensing  for  atmospheric  thermal-structure 
resolution  is  thus  based  on  processing  large  quantities  of  measurements  of 
channel  radiances.  The  CLRX  processing  objective  is  to  reduce  this  data  in 
volume,  distilling  out  the  Inherent  resolution  of  the  dear-column  radiance 
components  In  the  objective  scale.  The  diagnostic  model  assumes  that  each 
spot,  of  a  duster  of  spots  viewed  in  a  limited  cloudy  region ,  is  a  different 
mix  of  distinct  and  steady  radiance  components  ranging  anywhere  from  dear 
to  overcast  in  any  one  view.  "Distinct"  implies  that  the  model  compliance 
is  good.  "Steady"  implies  that  the  spatial  variabllitv  of  the  radiance 


components  is  small  compared  to  the  spot-to-spot  variabilities.  Not  much 
can  be  done  to  directly  enhance  the  desired  condition  of  distinctness.  But 
the  desired  condition  of  steadiness  can  be  enhanced  by  implicit  reduction 
of  all  variances  which  are  extraneous  to  the  cloud-mix  contrast  from  the 
measured  radiances  prior  to  model  diagnosis.  We  refer  to  these  operations 
which  are  intended  to  enhance  the  cloud- proportion  contrast  between  spots 
in  proximity  as  the  pre-processing  components  of  CLRBLCK. 

1.3.2  Limitations  of  the  SSE/DMSP  System 

Figure  1  is  an  example  of  the  scan  pattern  of  the  SSE/DMSP  system. 
Obviously  the  spots  are  much  too  sparse  to  warrant  diagnosis  of  cloud 
characteristics ,  local  to  each  spot,  by  Factor  Analysis .  The  spot  density 
is  not  adequate  for  anything  more  than  a  limited  two-state  interpretation 
of  the  cloud  structures.  Much  preprocessing  and  many  special  provisions 
are  required  to  distill  the  inherent  clear-column  resolution  from  the 
measured  radiances. 

Another  limitation  of  the  data  is  the  limited  significance  given  for 
each  measured  radiance.  The  significance  is  limited  by  the  number  of  bits 
in  transmission,  and  the  rounding  following  calibration  of  each  scanline. 

Such  a  system  characteristic  may  be  due  to  an  engineering  tendency  to  limit 
data  to  absolute  significance — not  appreciating  the  possible  values  of 
relative  information,  such  as  spot-to-spot  differences  in  the  present  context. 

The  diagnosis  components  of  models  CLRX  and  CLRXN  attempt  no 
more  than  a  two-state  interpretation  for  each  spot  of  the  scanline.  However, 
Imposition  of  a  two-state  interpretation  on  a  full  set  of  atmospheric  channels, 
even  though  the  adjustment-axis  components  are  individually  weighted,  is 
accompanied  by  characteristic  errors  when  the  cloud  structure  is  more 
complex.  On  the  average  these  errors  appear  to  negatively  bias  the  deeper 
components  of  the  diagnosed  clear-column  radiances. 


We  have  introduced  a  refinement  in  CLRBLOK  which  we  feel  should 
result  in  a  significant  increase  in  the  information  yield.  The  imposition 
of  the  two-state  interpretation  is  depth  limited,  to  the  optimum  validity  of 
the  interpretation.  The  diagnosis  makes  use  of  the  monotonic  character 
and  ordering  of  the  channels  according  to  the  depth  focus  of  their  atmospheric 
source  functions. 

Following  preprocessing,  channel  subsets — subset  1,2,3,  subset 
1,2, 3, 4,  subset  1,2, 3, 4, 5,  and  subset  1,2, 3, 4, 5,6 — are  each  assessed 
for  their  local  degree  of  compliance  with  a  two-state  interpretation.  At 
each  spot  the  optimum  subset  for  that  spot  is  selected.  The  components  of 
the  adjustment  axis  are  then  calculated  only  for  this  subset.  If  the 
optimum  atmospheric  subset  extends  to  include  channel  6 — implying  that 
the  whole  atmospheric  column  is  amenable  to  two-state  interpretation — 
then  adjunct  adjustment-axis  components  are  also  calculated  for  channels  7 
and  8  of  that  spot.  The  weight  for  each  calculated  axis  component  is 
calculated  according  to  the  local  variances  of  the  respective-channel 
component. 

Any  clear-column  radiance  components  produced  for  channels  7  and 
8  should  be  considered  to  be  byproducts  of  the  CLRBLOK  diagnosis.  Their 
reliability  is  subject  to  additional  considerations:  In  addition  to  being 
dependent  on  the  resolution  extending  to  channel  6,  the  accuracy  also 
depends  on  the  spot-to-spot  channel  uniformity  of  the  diagnosed  clear- 
column  radiances  for  each  of  these  channels:  Each  component  of  the 
adjustment  axis  is  based  on  the  assumption  that  local  spot-to-spot 
differences  are  due  only  to  differences  in  the  cloud  proportion  of  each 


1.4  The  CLRBLOK  Model 


1.4.1  The  Preprocessing  Components 

1.4. 1.1  Purposes 

The  primary  purpose  of  the  preprocessing  components  is  to  implicitly 
reduce  all  variations  inherent  in  the  measured  radiances  which  are 
extraneous  to  the  diagnostic  model — extraneous  to  the  assumption  that  local 
spot-to-spot  radiance  differences  are  largely  due  to  differences  in  cloud 
amount.  This  assumption  is  basic  to  the  assessment  of  clear-column 
radiances,  to  the  diagnosis  of  the  local  adjustment  axis,  and  to  imposition, 
in  the  total  information  blending,  of  a  degree  of  horizontal  continuity  for 
the  clear-column  atmospheric  radiance  components.  These  CLRBLOK 
refinements  extend  applicability  into  regions  where  surface-temperature 
variability  is  pronounced,  enabling  information  yield  for  the  atmospheric 
components  of  the  clear-column  radiances.  Additional  benefits  are 
indicated  in  Section  1.4. 1.4  below. 

1.4. 1.2  Local  Standard  Atmospheres  and  Transforms 

Nineteen  standard  pressure  levels  are  provided  for  specification  of 
first-guess  temperature  profiles:  0,  4,  7,  10,  20,  30,  50,  70,  100,  150, 
200,  250,  300,  400,  500,  700,  850,  925,  and  1000  mb.  The  temperatures 
specified  for  these  levels  are  subscripted  1  through  19  respectively: 
subscript  20  refers  to  the  surface  temperature.  The  temperature  profile  is 
specified  to  be  linear  in  p*  ,  where  k  =  2/7. 

First-guess  temperature  profiles  are  interpolated  from  available 
isobaric  temperature  fields  for  every  3rd  spot  of  each  scanline:  this  data 
includes  the  central  spot  numbered  13.  Temperatures  for  levels  not  now 
routinely  available  are  generated  by  linear-in-p*  for  T5  and  T~  , 
and  by  setting  T,  ,  T.  and  T  equal  to  T  ,  .  Sea-surface  temperatures  are 
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interpolated  from  available  analyses  for  all  spots.  These  interpolations 
involve  special  FNWC  subroutines  (see  Section  3.2). 

A  standard  atmosphere  is  defined  for  each  scanline — the  first-cuess 
temperature  profile  for  spot  13.  These  standard  profiles  are  generally 
continuous  from  scanline  to  scanline. 

The  available  transmission  functions  are  required  only  for  zenith-angle 
zero  degrees  and  zenith-angle  thirty  degrees .  These  transmission-function 
tables  are  expanded  by  interpolation  to  include  all  the  standard  pressure 
levels  listed  above.  Standard  radiances  are  calculated  for  each  channel  of 
each  scanline  for  each  of  the  two  zenith  angles. 

Based  on  a  linearization  of  the  Planck  function  about  the  star.card- 
atmosphere  temperatures  at  each  pressure  level,  transformations  are 
calculated  which  express  the  radiance  anomaly  for  any  channel,  v  ,  and 
zenith  ancle  0  or  30  degrees,  as  a  linear  combination  of  the  standard-profile 
temperature  anomalies,  including  the  surface-anomaly  contribution.  Several 
additional  parameters'are  calculated  for  each  channel  for  each  of  the  two 
zenith  angles:  a  mean  pressure  level  for  the  atmospheric  source  function, 
and  a  normalizing  factor  for  converting  each  radiance  anomaly  into  units  of 
temperature  anomaly.  The  relevant  formulations  for  these  transformations 
were  detailed  in  Holl  (1973a)  and  again  appeared  in  Holl  (1975a,  1975b  and 
1976a). 

Standard-atmosphere  values,  the  anomaly-transform  coefficients , 
atmospheric-source  pressure  levels,  and  temperature-normalization  factors 
are  calculated  only  for  the  zenith  angles  of  zero  degrees  and  thirty  degrees. 
Preprocessing  components  calculate  corresponding  values  ter  the  spot 
zenith  angles,  by  linear  interpolation  or  extrapolation  in  terms  of  an  optical 
path-length  equivalence. 


1.4. 1.3  Component  Operations 

All  measured  radiances  are  first  reduced  tc  anomalies  from  local 
s.ar,v.3.d  dea. -column  raoiances,  defined  for  each  scar.iine,  ir.cludinq 
zenith-angle  dependence.  Three  operations  which  act  implicitly0  on  the 
dear-column  radiance  components  follow:  Removal  of  the  surface  contribution 
from  all  measured  radiances  for  those  channels  which  penetrate  to  the  surface 
in  clear  conditions ,  by  proportional  subtraction  of  the  measured  window 
radiances,  normalization  of  the  radiance  anomalies  into  temperature— anomaly 
units,  ano  nadir  normalization.  Nadir  normalization  modifies  the  temperature- 
normalized  atmospheric  radiance  anomalies  by  interpolation  (downward 
extrapolation  in  the  case  of  cnar.r.el  6'  between  associated  source-function 
pressure-reference  levels  at  spot  zenith  angle,  to  the  corresponding 
pressure-reference  level  at  zero  zenith  ancle. 

These  operations  are  followed  by  subtraction  of  calculated  independent 
(i.e.,  first-guess)  temperature-normalized  nadir-normalized  radiance 
anomalies  from  standard. 

The  first-guess  radiances  are  calculated  using  only  the  required 
zero— zenith— angle  transtcrm3ticns .  The  independent  temperature-profile 
anomalies,  obtained  for  every  third  spot  of  each  scanlir.e,  are  transformed 
into  calculated  radiance  anomalies,  and  these  are  linearly  interpolated  for 
the  in-between  scots. 


o_  ,  ... 

-a oh  measurec  radiance  is  a  sum  of  clear-cclumn  and  cioud-to? 
components .  .An  implicit  operation  forms  linear  combinations  of  measured 
channel  radiances  all  of  one  spot,  in  order  to  effect  implicit  adjustments  of 
the  inherent  dear-column  components.  Corresponding  adjustments  :a!>;e 
place  in  the  cloud-top  components. 


Subtraction  of  these  calculated  first-guess  radiances  reduces  another 
of  the  variances  which  are  extraneous  to  the  diagnostic  model — the  spatial 
spot-to-spot  variability  of  the  inherent  clear-column  radiance  components. 
CLR3LOS  also  includes  provisions  to  exploit  these  calculated  independent 
radiance  estimates  in  the  diagnosis  components  of  Part  2. 

1.4. 1.4  Reversibility 

All  of  the  preprocessing  components  are  derived  from  transformations, 
developed  between  each  clear-column  radiance  and  the  temperature  profile , 
based  on  available  transmission  functions.  Accuracy  is,  of  course,  desirable 
to  make  the  preprocessing  components  as  effective  as  possible  in  the 
subsequent  diagnosis.  However  the  errors  introduced  thereby  are  not  of 
primary  concern:  The  preprocessing  components  can  be  reversed  after 
diagnosis  of  the  clear-column  radiance  components.  Reversals  are  performed 
in  the  restoration  components  of  CLRRLOK,  with  the  exception  of  two 
operations: 

(i)  The  diagnosed  clear-column  radiance  components  are  left 

nadir  normalized.  This  improves  utility.  Channel  calibration 
errors  have  only  a  second-order  effect  on  these  adjustments. 

(li)  The  surface  contributions  are  net  restored  to  the  atmospheric 
channels.  The  diagnosed  clear-column  window  radiances  are 
generally  not  as  reliable  as  the  diagnosed  and  weighted  clear- 
column  atmospheric-component  radiances.  The  diagnosed 
atmospheric  components  should  be  calibrated  and  exploited 
directly  for  use  in  atmospheric  thermal-structure  resolution 
enhancement. 


1.4. 1.5  The  Freprocessed  Radiances 


The  preprocessing  components  transform  the  measured  radiance  of 
channels  1  through  6,  into  a  corresponding  matrix  of  values,  in  several 
component  operations:  anomalies  from  local  standard  radiances,  implicit 
atmospheric  component  of  the  clear-column  component,  temperature 
normalized  equivalent  units,  implicitly  nadir-normalized  clear-column 
components,  anomalies  from  the  corresponding  independent  calculated 
(i.e.,  first-guess)  equivalent  quantities.  The  measured  channel  7  radiances 
are  transformed  into  temperature-normalized  anomalies  from  the  first-guess 
surface  temperatures .  The  measured  channel  8  radiances  are  unchanged  by 
the  preprocessing  components. 

This  three-dimensional  matrix  of  values ,  V  channels,  N  spots  and 
M  scanlines,  is  passed  onto  Part  2:  The  Diagnosis  Components .  The 
corresponding  independent  first-guess  estimates  of  the  clear-column 
components  are  zero  by  definition  for  channels  1  through  7. 

The  calculated  independent  (i.e.,  first-guess)  radiances  and  all 
ether  quantities  required  for  reverse  transformation  of  the  clear-column 
components  diagnosed  by  Fart  2  are  saved  for  use  in  Fart  3:  The  Restoration 
Components . 

No  matrix  of  associated  weights  is  attached  to  the  calculated  first- 
guess  radiance  estimates.  Such  associated  weights  are  dependent  on  the 
variances  inherent  in  the  analyzed  temperature  fields,  in  the  interpolation  to 
spot  locations  in  space  and  time,  and  in  the  transformation  coefficients  which 
in  turn  are  dependent  on  the  accuracy  of  the  available  transmission  functions. 
Instead,  when  exploited,  associated  weights  are  estimated  and  specified 

4 

only  as  function  of  channel  number  and  cycle  number. 


The  diagnosis  and  blending  components  are  repeated  three  times  ,, 
defining  three  cycles ,  in  order  to  purify  and  refine  the  weighting  of  the 
diagnosed  information  elements. 
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1.4.2  The  Information  Diagnosis  and  Blending  Components 


1.4.2. 1  Diagnosis  and  Blending  in  Terms  of  Adjustment  Values 

The  diagnosis  and  blending  components  are  formulated  in  terms  of 

adjustments,  o  ,  that  are  to  be  added  to  the  preprocessed  re  ances  , 

v  /  n  / 

E  _  ,  to  produce  the  equivalent  dear-column  radiance  components. 
v,n,^ 

The  subscripts  refer  to  channel  number,  spot  number,  and  scanline  number, 
respectively. 

As  is  the  practice  in  most  FIB  methodology  applications  ,  CLR3LCK 
includes  the  standard  optimum  number  of  three  cycles  through  the  information 
diagnosis  and  blending  components.  Three  cycles  allow  for  an  initial,  and 
two  progressively  more  stringent  evaluations  of  individual  information 
elements.  The  notation  of  denoting  the  resultant  values  of  the  blending  of 
all  assembled  information  by  a  superscript  asterisk  is  also  standard.  After 
each  cycle  of  diagnosis  and  blending  the  resultant  adjustment  values  are 
added  to  the  preprccessed  radiances  to  produce  the  resultant  clear-column 
radiance  equivalents: 

E*  m  =  E  „  +  O*  „  (I) 

v,n,^  v,n,^  v,n,”i 


To  begin  with,  the  only  estimates  that  are  available  for  ccmparitive 
evaluations  of  diagnosed  information  elements  are  the  independent  first-cuess 
estimates.  These  translate  into  the  adjustment  estimates 


O  =  -  E 

v.n,’’!  v,n,^! 


(2) 


But  each  subsequent  blending  operation  produces  new  estimates  —  mere 
reliable  and  hence  mere  effective  for  comparitive  evaluations  in  the  diagnosis 
component  of  the  following  cycle. 
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1.4. 2. 2  The  Diagnosed  Information  Elements 

Three  types  of  information  elements  are  diagnosed,  assembled,  and 
imposed  in  the  subsequent  blending: 

(i)  Absolute  Estimates . 

0 

Tentative  weights,  A  ,  are  assigned  to  those  radiances  which 

are  assessed  to  be  cloud  free  on  the  basis  of  local  variances  in  measured 

radiances.  The  associated  adjustment  estimate  is  o  =  0  . 

'  v,n,n 

These  information  elements  are  assembled  with  the  independent 

(c) 

first-guess  estimates,  -E  .  which  are  given  the  weights  A 

v,n,»J  *  v,n,^ 

The  superscript  c  in  parentheses  denotes  cycle  number.  The  resultants  of 
this  assembly  are  given  by: 


-  A(C)  E 

v  ,  n , m. _ v ,  n .  m 

A<C)  r  +  A* 
v,n  ,m  v,n,-n 


(3) 


A 


v,n,*i 


v ,  n,"i 


(4) 


0 

The  tentative  weights,  A  ,  are  reassessed  in  each  cycle, 

v  ,  n  , 

making  use  of  the  first-guess  information  in  cycle  one,  and  of  the  preceding 
blended  resultants  in  cycles  two  and  three. 


The  weights  given  to  the  first-guess  estimates  should  be  low  in 
cycle  one  because  of  calibration  uncertainties  and  because  the  estimates 
are  not  independent  from  channel-to-channel  and  spct-to-spot.  The  weights 
are  reduced  further  in  cycle  two,  and  are  made  negligible  in  cycle  three  in 
order  to  leave  the  final  resultants  free  of  any  direct  contributions  from  the 
first-guess  estimates.  The  first-guess  information  is  exploited  primarily 
for  and  in  the  diagnosis  of  information  elements  inherent  in  the  measured 


radiances,  as  for  example  in  situations  where  ambient  channel  variances 
are  low  not  because  the  region  is  cloud  tree  but  because  of  a  solid  uniform 
cloud  top. 

(ii)  Adjustment  Axes . 

A  weighted  set  of  difference  vectors  is  formed  for  each  spot  n,n  ,  by 
differencing  its  set  of  measured  radiances  with  those  of  six  ambient  spots  — 
four  on  the  same  scanline  and  the  corresponding  two  spots  from  adjoining 
scanlines.  The  best  available  adjustment  vector  for  the  spot  forms  a 
seventh  vector,  suitably  weighted  according  to  cycle  number. 

The  coherence  of  these  seven  vectors  is  evaluated  in  subsets  of 
channels :  1,2,3;  1,2, 3,4;  1,2, 3, 4, 5;  and  1,2, 3, 4, 5, 6.  The  subset  of 
best  coherence,  if  the  best  is  deemed  adequate,  is  then  combined  to  form 
a  single  adjustment  axis  of  that  many  channel  components.  Component 
weights  are  calculated,  based  on  variances  within  the  relevant  subset  and 
local  sample  of  difference  vectors.  If  the  optimum  subset  includes 
channel  6,  then  adjustment  components  are  also  formed  and  weighted  for 
channels  7  and  8. 

(iii)  Horizontal  Continuity 

The  spot-to-spot  continuity  in  the  clear-column  radiance  components 

implies  that  there  is  little  difference  in  the  resultant,  E  +  e* 

v ,  n ,  ~  v ,  n , "! 

between  spots  adjoining  in  a  scanline  and  adjoining  between  scanlines. 

The  corresponding  weights  assigned  to  these  general  information  statements 
are  specified  to  be  largest  in  cycle  one,  thus  forcing  the  first  blending 
resultant  to  produce  a  larger  objective  scale  of  resolution.  Although 
representative  of  a  larger  scale  than  the  final  objective,  the  first-cycle 
resultants  are  relatively  more  reliable  and  hence  are  useful  in  evaluatinc 
information  elements  in  cycle  two.  The  subject  continuity  weights  are 
reduced  in  the  cycle-two  blending,  and  again  in  cycle  three  to  allow  the 
desired  inherent  scale  of  resolution  to  emerge. 


1.4. 2. 3  The  Blending  Operation 


In  each  cycle  all  assembled  and  weighted  information  elements  are 

blended  to  produce  the  optimum  resultant  adjustment  field,  c  .  This 

v,n," 

resultant  minimizes  the  associated  error  functional  of  the  information  system. 
Minimization  amounts  to  solution  of  the  associated  system  of  simultaneous 
equations: 


MO*  =  F  (5) 

Sfc  -*•  — 

where  the  solution  vector  o  *  and  the  forcing  vector  F  have  an  element  for 
each  radiance.  The  matrix  M  is  symmetric,  positive  definite;  its  elements 
are  formed  entirely  from  information  weights. 

The  solution  is  developed  by  the  method  of  Blending  by  Weighted 
Spreading  (BWS)  shifting  to  Successive  Over  Relaxation  (SOR)  whenever 
and  wherever  the  blending  weight  approaches  one  for  any  element  of  the 
v.n.’J  array. 


1.4. 2. 4  The  Associated  Resolution  Weights 

Associated  resolution  weights,  A*  ,  are  required  for  the  final 

v,n,^ 

resultant  clear-column  estimates ,  O  ,  produced  bv  the  third  cvcle  of 

v,n,n 

blending.  According  to  the  FIB  methodology  the  associated  weights  appear, 
inverted,  as  the  diagonal  elements  of  the  inverse  of  matrix  M  of  Eq.  (5). 
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Very  little  continuity  coupling  is  imposed  between  scanlines  in  cycle 
three.  Hence  reasonable  approximations  of  reliability  can  be  calculated  for 
each  scanline  by  inverting  only  the  coefficient  matrix  for  that  scanline. 


1.4.3  The  Restoration  Components 


The  restoration  components  simply  reverse  corresponding  operations 
of  the  preprocessing,  with  the  exceptions  already  mentioned  in  Section 
1.4. 1.3.  The  reversals  are  performed  explicitly  in  terms  of  the  diagnosed 
clear-column  components . 

1 . 5  Thermal -Structure  Blending 


We  have  appended  a  Thermal-Structure-Blending  fTSB)  capability  to 
the  CLRBLOK  job  stream  in  order  to  present  a  more  complete  functional 
capability.  This  version  demonstrates  how  efficiently  the  relevant  information 
in  clear-column-radiance  estimates  can  be  used  to  directly  enhance  the 
resolution  of  the  atmospheric  thermal -structure. 

The  TSB  capability  is  less  comprehensive  than  is  the  mass -structure 
blending  scheme,  named  STRBLND,  mentioned  in  Section  1.2.1.  3TRBLND 
involves  more  than  the  thermal  profile;  included  are  provisions  for  height 
estimates  of  pressure  levels. 


It  should  be  noted  at  the  outset  of  this  discussion  that  there  is  a 
reversal  of  objectives  between  CLRBLOK  and  TSB.  Relationships  between 
channel  radiances  and  the  thermal  profile  are  exploited  in  both  schemes . 

In  the  case  of  CLRBLOK  the  relationships  are  used  to  exploit  information 
concerning  the  thermal-structure  for  its  relevance  to  the  clear-column 
radiance  resolutions,  and  in  the  case  of  TSB  the  relationships  are  used  to 
exploit  information  concerning  the  clear-column  radiances  for  its  relevance 
to  the  thermal- structure  resolution. 


In  CLRBLOK  the  objective  is  maximum  resolution  of  the  clear-column 
radiances.  In  TSB  the  objective  is  maximum  resolution  of  every  element  of 
the  thermal-structure  model.  [The  thermal-structure  model  is  that  described 
in  paragraph  I  of  Section  1.4. 1.2.  ]  The  same  linearized  approximations  are 
used  to  express  relevant  independent  (i.e.,  first-guess)  estimates.^  Beth 
applications  of  information  blending  must  be  carried  out  to  effect  both 
objectives.  The  end  results  are  not  equivalent.  Appreciation  of  the 
difference  is  essential  to  the  understanding  of  how  the  FIB  methodology 
derives  maximum  relevance  from  available  information. 

Relating  a  clear-column  radiance-anomaly  estimate  to  a  linear 
combination  of  the  temperature-anomaly  elements  of  the  thermal-structure 
profile  model  involves  two  contributing  error  variances:  the  error  variance 
associated  with  the  diagnosed  clear-column  radiance,  and  the  error  variance 
inherent  in  the  form  of  the  linearized  anomaly  relationship. 

In  order  to  demonstrate  the  TSB  capability  we  have  used  the  sets  of 
linear-transform  combination  coefficients  which  are  computed  by  preprocessing 
component  Cl.l  for  scanline  5  of  the  development  sample  (see  Sections  2.1 
and  2.2).  The  diagnosed  clear-column  radiance  weights  are  reduced  rather 
sharply  by  addition  of  large  estimates  of  the  variances  attributed  to  the 
transforms . 

g 

A  subroutine  is  required  for  introduction  into  the  CLRBLOK  job  stream. 
Its  function  would  be  to  replace  the  available  standard  coefficients  by 


wIn  CLRBLOK  we  choose  to  minimize  the  direct  contribution  of  the 
calculated  clear-column  radiances  in  the  diagnosed  resultant  clear-column 
radiances,  in  order  to  use  these  diagnosed  radiances  as  independent 
satellite  information  for  subsequent  enhancement  of  the  thermal-structure 
profile. 

^This  requirement  is  included  among  our  recommendations  (see 
Section  4.) 
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coefficients  which  are  optimized  for  the  local  conditions.  These  transformations 
should  be  statistically  refined  in  order  to  minimize  the  associated  error 
variances . 


The  thermal-structure-blending  capability  includes  the  following 
features: 

The  blending  is  performed  in  terms  of  potential-temperature  values 
for  the  standard  pressure  levels .  The  independent  temperature 
estimates,  T.  ,  and  associated  weights,  A.  ,  are  transformed 
accordingly. 


The  clear-column  radiance  estimates  are  transformed  into  estimates 
of  linear  combinations  of  the  potential-temperature  values.  The 
weights  associated  with  these  estimates  are  reduced  in  accordance 
with  contributing  error  variances  ascribed  to  the  relationship  for 
each  of  the  six  atmospheric  channels. 

The  blended  solution,  9*  ,  and  associated  resolution  weights ,  A*, 
are  calculated  by  explicit  matrix  inversion. 

The  solution  profile  is  tested  to  be  statically  stable:  It  is  required 
that  the  potential-temperature  values  decrease  from  one  standard 
pressure  level  to  the  next,  beginning  at  the  top  of  the  atmospheric 
column.  If  any  level  fails  the  static-stability  test,  because  the 
potential-temperature  value  is  found  to  be  greater  than  that  for  the 
level  just  above,  then  the  value  is  replaced  by  that  for  the  level 
just  above.  This  procedure  is  followed  sequentially  from  top  to 
bottom.  The  adjusted  solution  is  then  optimized,  within  the  constraint 
of  positive  static  stability,  by  sequential  SOR  passes  from  top  to 
bottom.  This  procedure  exploits  the  additional  information  that  the 
atmosphere  is  statically  stable  in  the  object  scale  of  resolution. 


The  potential-temperature  values  of  the  solution  are  transformed 
into  temperature  values,  T*.  The  associated  resolution  weights  are 
transformed  accordingly. 

Independent,  first-guess,  temperature  profiles  are  obtained  for  use  in 

CLRBLOK  for  every  third  spot  of  each  scanline.  The  Thermal-Structure 

Blending  (TSB)  capability  is  programmed  to  process  this  subset  of  spots 

7 

excluding  spots  1  and  25  of  each  scanline.  The  Sample  run  (see  Section  2), 
which  includes  TSB  applications  for  63  profiles  has  demonstrated  the  efficiency 
and  effectiveness  of  this  capability. 


Clear-column  radiances  are  not  diagnosed  for  spots  numbered  1,2, 
24  and  25  of  each  scanline  because  of  the  sparsity  of  spots  in  scanline 
extremes,  and  for  convenience  in  avoiding  boundary  complications  in  the 
programming . 
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2. 


DEVELOPMENT  OF  THE  JOB  STREAM 


2.1 

Part  1 


Part  2 


Component  Operations  of  the  Job  Stream 

The  Preprocessing  Components 
C1.0  Initialization 

Cl.l  Generation  of  Standard  Profile  Transforms 

Cl. 2  Local  Standard  Radiances,  and  Anomalies  from  Standard 

Cl.  3  Implicit  Removal  of  Surface  Contribution 

Cl.  4  Implicit  Temperature  Normalization 

Cl. 5  Implicit  Nadir  Normalization 

Cl.  6  First-Guess  Radiance 

Cl.  7  Optional  Linear  Transformation  of  Atmospheric  Channels 
(not  presently  used) 

Cl. 3  Subtraction  of  First-Guess 

Cl.  9  Organization  of  Scanline  into  Block  Format 

Cl.  10  Save  Labeled  Common  31ocks  for  Later  Use 

The  Diagnostic  Components 
C2.0  Initialization 

C2 . 1  Ambient  Variance  in  Preprocessed  Radiances 
C2.2  Adjustment  Axis  Estimates  for  Each  Spot 
C2.2.1  Initialization 
C2.2.2  Local  Contrast  Vectors 
C2.2.3  Optimum  Depth  of  Model  Endorsement 
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C2 . 2 . 4  Local  Adjustment  Axis 

C2 . 2 . 5  Adjustment  Axis  Weights 

C2.2.6  Selective  Print-Out 

C2.3  The  A  Weights,  and  Assembly 

C2.3.1  Initialization 

C2.3.2  Clear  Measure  Assessment 

C2.3.3  Assembly  of  Estimates  and  Weights 

C2.4  Blending  for  Resultant  Estimates 

C2. 4. 1-2. 4. 2  Initialization 

C2.4.3  The  31ending  Operation 

C2. 4. 4  Process  the  Solution: 

Cycle  1  and  2:  return  to  Section  2.2 
Cycle  3:  proceed  to  Section  2.3 

C2.5  The  A*  Approximations  for  each  Scanline 

C2.5.1  Fill  the  Inversion  Matrix 

C2 . 5 . 2  Invert  the  Matrix 

C2.5.3  A*  Values 

Part  3:  Restoration  Components 
C3. 0-3.1  Initialization 
C3.2  Add  First-Guess  Back 
C3.3  Re-Scale  Radiance  Anomalies 

C3.4  Add  Surface  Contribution  Back  (Optional — not  exercised) 

C3.5  Add  Local  Standard  Back 

C3.6  Re-Calibrata  Resolution  Weights 
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Part  4:  Thermal-Structure  Blending  Components 


C4.1  Initialization 

C4.2  Refinements  of  Coefficients  and  Specifications 
(Provision  for  future  subroutine) 

C4.3  Estimate  and  Weight  Transformation 

C4.4  Explicit  Blending  Solution 

C4.4.1  Fill  Blending  Matrix 

C4.4.2  Invert  Matrix 

C4.4.3  Form  Forcing  Vector 

C4.4.4  Resultant  Profile  9* 

C4.4.5  Associated  Resolution  Weights 

C4.5  Stabilize  the  Profile 

C4.5. 1  Test  and  Adjust  Profile 

C4.5.2  Optimize  Adjustments 

C4.5  Organize  Resultant ,  E** 


2 . 2  Print  Potions 

The  CLR3LOK  job  stream  has  been  designed  to  include  a  '.vide  choice 
of  printout  options  to  aid  in  the  initial  program  debugging,  in  developing 
familiarization  with  component  functions,  and  for  evaluating  program 
modifications  and  tuning  adjustments. 
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The  amount  of  printed  output  can  be  varied  for  each  scan  and/or 
block  of  scans: 

All  Sections: 

Numerical  values  for  each  scan,  and  the  block,  are  stored 
in  array  MOPT(20),  located  in  Labeled  Common  Block 
PARBL.  These  values  range  from  1  (for  complete  print-out) 
to  5  3  (no  print -out) . 

MO PT 2  .  . .  MOPTg  contain  print  options  for  scan  41. .  .49. 
MOPT1Q  contains  print  option  for  the  whole  block. 

MO  is  generally  =  the  current  print  option. 


Sections  4.3. .  .4.6: 

If  MO  =  3,  then  any  profile  spot  can  be  selected  for 
detailed  print-cut. 

Print  option  values  are  stored  in  arrays  MON? (7),  MOM(9), 
and  MONPM(7,3)  in  Section  4.0. 


These  details  are  covered  in  C4.  1  (Initialization)  of  the 
Program  Documentation. 


The  lists  of  Common  Blocks,  arrays,  and  other  values  with  their  accompanyin 
print  options  are  included  in  the  functional  description. 


V 


» 


I 
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2.3 


he  Development  Samnle 


The  sample  block  of  nine  scanlines  shown  in  Figure  1  was  used  as 
the  development  sample  for  program  CLR3LOK.  This  sample  block  was 
selected  rrom  the  available  collection  of  radiances  because  it  was  the 
longes-  continuous  sequence  of  scanlines  over  open  ocean.  The  block 
was  not  selected  with  any  particular  expectations  as  to  yield,  nor  were 
•he  cloud  characteristics  considered  in  the  select’oo. 

The  information  yield  which  is  expressed  by  the  reliability  weights 
of  the  diagnosed  dear-column  radiances  is  dependent  on  the  characteristic 
o:  dee  cloudiness.  A  clear  region  generally  gives  the  highest  yield,  ^he 
CLRSLOK  capability  is  designed  to  diagnose  the  inherent  information  in  the 
measured  radiances  no  matter  how  great  or  how  meager  that  yield  might  be. 
It  is  even  mere  important  that  raise  information  net  be  generated  and 
assimilated  than  it  is  that  good  information  be  produced. 

2 . 4  Preliminary  Ccmocnents 


The  input  requirements  for  Program  CLR3LCK  include: 

■a)  The  measured  radiances  organised  by  channels,  spots  and  scanlines, 
with  spot  locations,  zenith  angles  and  timing  data. 

vb,1  nSt-.ma.es  ot  the  temperature  profiles  for  everv  third  spo*  of  gagn 
iach  prorile  is  expressed  by  the  temperature  (”R)  a*  oac'"'  o*  1  3 
pressure  levels. 


(c)  3 


stimates  of  the  sea-surface  temperature  for  each  spot,  and  a  land/ sea 
.ta-:unct:cn  (i.e. ,  0  or  1)  value  for  each  sect. 


(d)  T 


transmission  runcticns  for  each  channel  for  zenith  ancles  0°  and  31'  . 
These  functions  are  used  in  Part  1  (subroutine  RADCCF)  for  theoretical 
calculation  cf  the  transforms  used  in  ,■»». 


radiances ,  and  are  used  in  Part  4  in  thermal-structure  b leucine .  'Me 
consider  this  dependency  on  theoretical  transmission  functions  to  be 
an  interim  necessity,  awaiting  replacement  by  superior  means. 

In  an  operational  context  the  required  input  would  pass  through  the 
processing  center  in  a  lagged  real-time  mode.  Per  present  development 
purposes  the  above  data  was  assembled  by  R&D  procedures.  These 
procedures  have  been  designated  Part  0  of  Program  CLRBLOK,  but  they 
should  not  be  considered  part  of  the  job  stream.  The  input  preparations 
for  CLRBLOK,  in  an  operational  context,  will  depend  on  the  specifics  of 
the  operatinc  system. 


2.5  The  Samole  Performance 


The  yield  derived  from  the  sample  block  of  scanlir.es  turned  out  to 
be  high  for  channels  1,  2  and  3,  effectively  nil  for  channels  5,  7  and  8, 
and  intermediate,  and  spotty,  for  channels  4  and  5,  all  accordinc  to  the 
product  associated  reliability  weights.  This  untuned  diagnosis  by  CLRBLOK 
implies  that  the  region  has  a  complex  cloud  structure  above  a  low  overcast. 
The  comolex  cloudiness  has  been  confirmed  bv  examination  o*  *"~e  co’'ourr,ar* 
synoptic  data.  Figure  2  is  a  surface  analysis  with  station-model  plots. 

The  four  stars  in  Fig.  2  indicate  the  corner  points  of  the  cuasi-rectar.cular 
region  covered  by  the  nine  scanlines.  Considering  the  synoptic  situation, 
the  atmospheric  depth  or'  model  relevance  and  the  resuiting  information 
yield  of  the  CLRBLOK  diagnosis  are  understandable .  The  medal  performed 
according  to  performance  expectations  . 
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2.6 


Proem  re  Tlnin  or- 


The  following  timing  figures  were  determined  far  the  CLRSLOK  ich 
stream,  including  the  thermal-structure-blending  supplement  in  CLRBL3. 
These  figures,  in  C?  seconds  of  execution  time,  apply  to  the  processing  of 
the  sample  block  cf  nine  scanlines,  using  FNWC's  R&D  CDC-5500  computer 
named  "HAL". 


With  Printout 
(as  delivered) 

Without  Printout 
but  including  all 
print  decision 
steos 

CLRBL1 

16.541 

2 . 190 

CLRBL2 

1 33.913 

158.652 

CLRBL3 

22 .999 

17.091 

Total 

223.453 

137.933 

Removal  of  all  print  provisions  would  make  a  s  i  rmticant  roducticr 
in  the  length  of  the  program  listing--about  a  30k  reduction. 


2 . 7  Program  Documentation 


An  extensive  program  documentation  of  the  CLR3LO\  :ob  stream, 
including  the  preliminary  components  and  the  thermal-structure-olendmc 
capability  supplement,  has  been  prepared  and  assembled  under  separate 
cover.  This  volume  of  documentation  includes  th  ■*  rollowing  material: 

General  flow  diagrams  showing  structure  e:  complete  'Ob  stream, 
processing  programs,  and  accompanying  subroutines. 


Detailed  functional  description  of  complete  job  stream: 
Preliminary  Components  ,  Preprocessing  Components,  Diagnosis 
Components,  and  Restoration  and  Thermal-Structure  Blending 
Components . 

Control  card  listings  for  each  job-stream  unit  are  included  in  the 
functional  description,  as  well  as  in  each  file  listing. 

Detailed  common  block  definitions. 

Subroutine  calls,  and  parameters  used. 

Print  Options . 


Program  listings 
and  execution  printouts) 


and  the  sample  production  run  vincludinc  compilation 
have  also  been  delivered  to  the  sponsoring  acency 
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FORMULATIONS 


3.  1  Ifovelormor.ts  ,>»'  ;h.-.  ~y  >r.d'.r.q  S  vs _m 

3.1.1  The  Functional  and  the  System  of  Blending  liquations 

Applications  of  the  Fields  by  Information  Blending  (FIB1  methodology 
begin  with  the  design  of  the  functicnal--with  the  design  of  provis:  -ras  in 
the  functional  to  include  the  various  available  relevant  terms  or  information. 
The  functional  includes  various  elements  of  information,  and  these  elements 
with  their  associated  weights  serve  as  repositories  for  the  assembly  of 
information  gathered  in  forms  of  weighted  estimates.  The  CLRX,  CLRXX, 
and  CLRBLOK  functionals  include  a  form  of  information  which  is  unique  in 
FIB  applications:  the  components  of  adjustment  ox  v-  for  subset  vectors  of 
the  object  array.  The  functional  expresses  the  sum  of  the  weighted  squared- 
disparities  to  all  assembled  elements: 
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Notation: 


rne  subscripts  v ,  n , ' 


**  &  -  *  ■  r  -  y 


the  channel,  spot,  and  scanlir.e  number 
respectively*  The  summation  sxtends  >v«r  the  entire  block  •  scanlinei 


I 

L- 


_  .  is  the  object  adjustment  that  Is  t  •  be  •  Ided  to  the  pre processed 
radiance,  E  „  ,  to  produce  the  equivalent  dear-column 

\  t  Vi  , 

radiance  components,  e*  -  E  _  .  for  channel  v  or  spot 

*  *  **  /  * 

of  scan  line  “  . 


c,.  _  „  is  the  assembled  direct  estimate  of  o*  and  A 
^  -  n  <  •  v  ,  n ,  n  v ,  n , ' 

weight . 


3  is  the  weight  of  the  estimate  that  there  is  cere  difference  between 

the  clear  components  of  spots  v,n,~  and  v,n-l,~  . 


w  ‘ s  the  weight  ot  tne  estimate  thet  tnere  n c t* c  ciizfv-r  ^  •-'n 

,  n  /  ♦ 

the  clear  components  of  spots  v,n,n  and  v , n , ~-i  , 


„  is  the  v  component  of  the  adjustment-axis  estimate  at  the  o,~ 

v  » n » 

spot. 

_  is  the  unknown  length  factor  d  the  adjustment-axis  estimate  at  th 
n,"  spot. 


'•V  _  is  the  weight  ccrresoondtnc  to  the  L 

v  f  n  f  • 


\  _  sdfustnont  -r.at?, 

t*1 


*  . 


The  equation  for  S  ,  Eg .  (11),  car.  be  considered  to  be  inhere".: 

v ,  n ,  ~ 

in  Eq.  (10).  Equation  (10)  should,  first  of  all,  be  recognized  as  a  weighted 

summation  of  several  estimates  ’■v'°  quanritv  c  ^rd  3 

v ,  n ,  ~  '  v ,  n ,  n 

should  be  recognized  as  the  sum  of  the  weights.  The  individual  estimates 
appear  to  the  right  of  each  heavy  dot,  and  the  weights  are  the  coefficient 
factors  to  the  left  of  each  heavy  dot.  We  propose  that  weighted  combination 
of  estimates  be  written  using  the  notation  of  heavy  dots  to  separate  individu; 
estimates  and  their  coefficient  factors.  And  we  further  propose  that  a  socor.c 
equation  is  implicitly  specified  by  such  no  tatter.-- -the  equation  obtained  by 


dropping  each  estimate  to  the  right  of  each  dot,  summing  only  the  coeffioiei 
weights.  According  to  this  notation  and  convention,  Eg.  i.11)  is  implicit  in 
Eq.  (10). 
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